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We solve the standard Lifshitz-Slyozov (LS) model with conservation of total particles in the limit 
of small super-saturation. The new element is an effective initial condition that follows from the 
initial exhaustion of nucleation as described in a previous paper The effective initial condition 
is characterized by a narrow distribution of cluster-sizes, all much larger than critical. In the 
subsequent solution, one of the LS similarity solutions emerges as the long-time limit, as expected. 
But our solution tells more. In particular, there is a "growth" era prior to what is usually called 
"coarsening." During "growth" the clusters (all of nearly the same size much larger than critical) 
eventually exhaust the super-saturation (the exhaustion of nucleation in the previous era results 
from only a small decrease in super-saturation). This allows the critical size to catch up to the 
clusters, and the traditional "coarsening" begins: Subcritical clusters dissolve and fuel the growth 
of the remaining super-critical clusters. Our analysis tracks the evolution of cluster sizes during 
growth and coarsening by complimentary use of asymptotic and numerical methods. We establish 
characteristic times and cluster sizes associated with growth and coarsening from physical parameters 
and the initial super-saturation. The emerging distribution is discontinuous at the largest cluster 
size, and thus selects the discontinuous LS similarity solution as the long-time limit. There are 
strong indications that the smooth similarity solution proposed in the original LS paper emerges on 
a, yet longer, "late-coarsening" time-scale. 

PACS numbers: 81.10.-h, 68.43.Jk 



Introduction 

We analyze the growth of clusters in a monomer bath 
with conserved total monomer density in the small super- 
saturation limit. The description of late stage coarsening 
according to the classic Lifshitz-Slyozov (LS) theory is 
well known in the literature of aggregation: The number 
of monomers in the largest clusters increases linearly in 
time, and the density of clusters shrinks as the smaller 
clusters dissolve back into monomers. But how docs it all 
begin, starting from pure monomer at "time zero"? In a 
previous paper [l[ the authors predict the cluster size dis- 
tribution that emerges from the nucleation process, and 
its long-term limit. In the current paper, the solution of 
the LS equations, with that long-term limit used as an 
effective initial condition, is tracked all the way to late 
stage coarsening. In this way, we obtain a "big picture" 
of the whole aggregation process, based on classical mod- 
eling ideas due to Becker-Doring 0], Zeldovich and 
Lifshitz-Slyozov 

Here is the summary using Fig. ([1]) as a visual guide. 
The horizontal axis is the largest cluster size, n,„, the 
vertical is time t, both with logarithmic scales. This 
graph of t vs. rim is based on the quantitative solution 
of the complete model. The plane is divided into hori- 
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zontal time-slices, "Creation" , "Growth" , and "Coarsen- 
ing." The characteristic time \t\ to exhaust nucleation 
is exponentially large in the initial free energy barrier 
G* against nucleation, with \t\ oc exp (iG^/fcsT). The 
timescale [t] is the thickness of the creation time-slice. 
In this time, the supersaturation undergoes only a small 
relative decrease and the initial clusters continue rapid 
growth. For diffusion limited growth in (nearly) con- 
stant super-saturation, the number of monomers n{t) in 
a cluster grows at a rate proportional to cluster radius, 
so ri cx m}^^ and it follows that n{t) cx t^^"^. In particular, 
the asymptotic line in the creation time-slice of Fig.[l]has 

2/3 f 

a 2 : 3 slope consistent with n„i oc -^j. In this way we 
see that the characteristic size [n] of clusters during the 
creation era is proportional to [t]'^/^ cx exp {^G^^/kBTy 
The width of the cluster size distribution grows more 

1/2 

slowly, like {t/[t]) ' . Hence, the relative width of the 
cluster size distribution becomes small during the tail of 
the creation era, oo. The actual profile of this nar- 

row distribution is determined from the time history of 
the nucleation rate per unit volume j{t), derived in 

In the next time-slice, labeled "Growth" , the nearly 
homogeneous population of rapidly expanding clusters 
seriously depletes the super-saturation. This depletion 
causes their rapid growth to stop when their (common) 

cluster size reaches -jr- Here, e with < £ <C 1 is the 
initial super-saturation. In Fig.[l] this arrested growth is 
represented by the nearly vertical segment. The charac- 
teristic time which measures the thickness of the growth 
time-slice is 
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FIG. 1: The graph of time vs. maximal cluster-size in logarithmic scale. The graph 
shows two obvious regimes (creation and coarsening) separated by a "kink" in the 
graph (growth and "pre-coarsening" ) The scales n,, [n] and (T^[n]/_Re^ of cluster size 
are marked in the plot. 



The next era, coarsening, begins when the critical clus- 
ter size, much smaller than the characteristic size during 
creation and growth, "catches up" with the clusters' size. 
Clusters smaller than critical shrink, and the monomers 
they shed are taken up by the larger, growing clusters. 
Thus, the distribution widens. The characteristic size of 
the clusters during the coarsening era remains the same 
as it was during the growth era, but the timescale is much 
longer: [n]3[<]/e^. 

In late stage coarsening, t ^ [7i]3[t]/e2^ ^j^^, dygter 
size distribution asymptotes to a similarity solution of 
the LS equations, and we finally reach the stage when 
rim is linear in time. Indeed, the asymptotic line in the 
coarsening time-slice of Fig. [T] has the characteristic 1 : 1 
slope. 

Thus concludes the "brief history" of aggregation ac- 
cording to the classical ideas of BD, Zcldovich, and LS. 
We highlight some collateral results. First regarding time 
and size scales: By introducing a physical initial condi- 
tion that represents the initial nuclcation process, we ul- 
timately determine the characteristic time to reach coars- 
ening and the characteristic cluster size, as functions of 
the physical parameters and initial supersaturation. In 
particular the time to reach coarsening, [t]c = [n] 3 
is exponentially large in the initial free energy barrier 
G*, even relative to the time [t] of the creation era. An 
actual physical time-scale cannot result from the LS the- 



ory alone. Indeed, the famous similarity solutions result 
precisely from scale invariance. 

This brings us to a peculiar detail: The late-stage 
coarsening similarity solution that is selected by our so- 
lution of the LS equations is discontinuous at the largest 
cluster size. It is widely believed that the physically cor- 
rect similarity solution is the smooth, C°°, one. In the 
discussion section we propose that during an additional 
era following coarsening, the distribution evolves further 
and tends to the smooth C°° similarity solution. 

The organization of the paper is as follows: In Sec- 
tion [H we review the origins of the cluster size distribu- 
tion during the tail of the creation era, as set forth in 
This, of course, is the effective initial condition for 
the growth era, treated in Section |TT1 Section HITl treats 
the coarsening era and its asymptotic matching with the 
tail of the growth era. Here there is an additional twist: 
The coarsening era is evolved numerically, whereas the 
effective initial condition inherited from the growth era 
comes from an analytic solution. The switching from 
analytic to numerical solution is controlled as a func- 
tion of the numerical resolution, so we have a de-facto 
"analytic-numerical matching." One corollary of this ex- 
panded sense of matching is the analytic determination 
of a time delay for the onset of coarsening, proportional 
to logi. 

In Section [IVI we re-derive the family of similarity so- 
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lutions using our notation. The discontinuous member 
of this family is determined and the convergence of the 
numerical solution to it is verified. 



I. THE PHYSICAL MODEL AND EFFECTIVE 
INITIAL CONDITIONS 

In the classic Lifshitz-Slyozov (LS) theory, the number 
of monomers n = nit) in a cluster satisfies the ODE of 
diffusion limited growth: 



- cr), (3(47r)2) Dv^^ fs 



(1.1) 



Here, ?/ is the chemical potential of monomers in the bath 
(in units of kgT) relative to monomers in the bulk of 
clusters. When the monomer density /i approaches the 
saturation density fa, for which the monomer bath would 
be in equilibrium with an "infinite" cluster, we have the 
asymptotically linear relation 



_ h-fs 



(1.2) 



In (|l.ip . a is the dimensionless surface tension constant 
so that the interfacial free energy of a cluster with n 
monomers is ^n^aksT. In the definition of the rate 
constant V, D denotes the diffusivity of monomers in 
the bath and v is the monomer volume inside clusters. 

The state variable of the LS equations is the cluster- 
size distribution r(n, t), so that the density of clusters 
with size n between ni and 712 is J^^ r{n, t) dn. There 
are two basic equations: First, the convection PDE 



dtr + dn \t> (rin3 - a'j rj 



= 0, for 71 > 0, 



(1.3) 



represents transport of clusters in the space of their size 
n by the diffusion limited growth "velocity" in (jl.ip . Sec- 
ond, the conservation of monomers couples the value of 
the super-saturation, 7/ and the solution. The conserva- 
tion of monomer is expressed approximately by 



/ = (1 + 7/)/s + / nr{n, t) dn 



(1.4) 



Here, the total monomer density /, a constant in time, is 
the sum of monomer density /i = (1 + 77) (from (|1.2p ) 
in the bath, and the the density of monomers in clusters 
is approximated by the integral. 

In the convection PDE p.3|l . a is positive, so charac- 
teristics in the (n,t) plane are absorbed by the 7i— axis. 
Hence, the 71— axis is a "sink" , representing the com- 
plete dissolution of subcritical clusters. This is consis- 
tent with the assumption that creation of new clusters 
by fluctuation over the critical size is negligible during 
the "growth" and "coarsening" eras. In a previous pa- 
per [l| we derive scaling units [t] , [r] , [77] of time t, cluster 
size 71 and cluster size density r that characterize the cre- 
ation era. It is convenient to express the characteristic 
scales of the growth and coarsening eras as multiples of 



these creation era scales. Hence, we carry out a prelimi- 
nary non-dimensionalization of (|1.31 II. 4p based on [t] , [r] , 
and [71]. The unit of chemical potential 77 is [7;] = e, the 
initial value of chemical potential in the pure monomer 
bath, before nucleation. The dimensionless equations are 



dtr + dn |7i3r — s| = 0, 



7/ = 1 



I r dn. 



in 71 > 0, (1.5) 
(1.6) 



In (|1.5[) . s is the scaled surface tension, exponentially 
small as e — ^ defined in appendix [Al in (|A.4p . Equations 
(|1.5| II. 6p are solved for r(77, t) subject to an effective ini- 
tial condition that arises from asymptotic matching with 
the creation era. In the previous paper we showed that 
at a range of time t, after nucleation is exhausted, but 
before the effects of growth change the super-saturation 
significantly, r{n, t) is asymptotic to a narrow distribu- 
tion is approximated by 



1 ■ ( N-n 
r{n, t) = <! N^/^ 
0, 



< N -n< 



otherwise. 

(1.7) 

Here, 71 = N{t) is the size of the largest cluster, approx- 
imated by 



N{t) 



^<)^ (fort = 0(1)). 



(1.8) 



The function j{t) is the dimensionless nucleation rate 
whose graph is shown in Fig. [21 In our previous paper 
Q it is shown that j{t) satisfies the integral equation 

l0gj(t)=- f\l{t-T))\j{T)dT. (1.9) 



Its solution, j{t), decays to zero faster than exponential 
as t 00. For t > 5, j{t) is a. negligible fraction of its 
initial value, so we truncate the support ofj{t) toO < t < 
5 for simplicity of presentation. This explains the upper 
limit 5N3 in (|1.7p . The reduction in monomer density 
due to aggregation during the creation era is denoted R 
and is calculated from j: 



R 



j{T)dT w 1.7117. 



(1.10) 



The value R ~ 1.7117 is, of course, a scaled density. To 
get a physical density one needs to multiply it by [77][r], 
with [71] and [r] given by (|A.2ilA.3p . 



II. GROWTH ERA 

During the growth era, the cluster distribution is still 
approximated by (|1.7p . but the growth of the largest 
cluster-size N{t) slows relative to the t^^"^ growth law 
due to the depletion of supersaturation. Here is a 
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FIG. 2: The solution of Eq. (|1.4|l is the rate of production of 
new clusters j = e''''''' during the nucleation era. The tail of 
the flux rate decays super-exponentially in t. The dashed line 
is the asymptotic result found in [J]. 



brief summary of the argument. In the convection PDE 
(|T5p . the component rj-n}/^ of convection velocity is much 
greater than one, so the scaled surface tension s is asymp- 
totically negligible. The convection PDE thus reduces 
asymptotically to 



dtr + d, 



= 0, in 71 > 0. 



(2.1) 



The corresponding physical idea is that most of the clus- 
ters are much larger than critical. It follows from (|2.ip 
that n^/^r{n, t) is constant along characteristics that sat- 
isfy 



(2.2) 



In (|2.2|1 . T] = ?7(t) decreases from (near) 1 in the begin- 
ning of the growth era to (near) at the end in a manner 
consistent with the conservation identity (|1.6p . We see 
that the characteristics determined by \2.2\ are contin- 
uations of the creation era characteristics, carrying the 
same values of n^^^r. 

This indicates a very simple construction of the asymp- 
totic solution for r(n, t) during the growth era. The de- 
tails are in Appendix |B] In summary, r(n, t) is concen- 
trated in a narrow front near the largest cluster size N , 
and there approximation (|1.7p applies. What changes is 
the evolution of iV(i), now described by the ODE 



iV = iV3 1 



N 
iVo 



4- (^-^^ 



Here, 1 — is the value of rj[t) consistent with the con- 



servation identity pTel) . The solution to ODE with 
N{Q) = and N{t) > for i > is given imphcitly by 



2 

log 



N 
iVo 



(2.4) 



where are the cube roots of 



-l,r2 



1: ro = e'3,ri = 
Figure [3] shows this solution as a "world 




FIG. 3: The "world-lines" of clusters created at the origin. 
The density of the lines corresponds to the density of clusters 
at each point. The length of the horizontal line from {N{t), t) 
to {No, t) (in units of No) is the supersaturation rj (in units 
of £.) 



line" in the (n, t) plane (dark line). The shaded area rep- 
resents the front where r(n, t) is concentrated. In the 



limit 1 < t < TV, 
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reduces to t 



3 7^2/3 
2 5 

from the creation era. 

2/3 



agreement with results 

In the opposite limit t ^ Nq''^, which corresponds to 
the tail-end of the growth era, N(t) asymptotes (expo- 
nentially) to the constant value A^o = /e^R. The size 
distribution asymptotes to 



A^ 



(2.5) 



independent of time. It is still narrow, with its support 
is concentrated in an interval of n with < Nq — n = 

0(7Vo*)« ATq. 

Why does the size distribution "stop dead in its 
tracks" ? In Fig. [Sj the length of the horizontal line seg- 
ment from {N{t)/NQ, t) to (1, t) represents the super- 
saturation 7] (in units of e) at time t. It asymptotes to 

2/3 

zero for t ^ A^q j ^^^^ the truncated convection veloc- 
ity rjn^/^d vanishes with it. The clusters "use up" the 
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super-saturation that fuels their growth. 

In summary, during the growth era, the clusters grow 
in a relatively narrow distribution until they reach a max- 
imal cluster size n = No (in units of [n]). The width of 

1/3 

the distribution is proportional to A^q . The time-scale 
of the era is A^q^'^ (in units of [t]), and roughly 10 of these 
time-units are needed for the narrow, stationary distri- 
bution to be established, as seen in Fig. [3l The growth 
of the clusters is fueled by the supersaturation, which 
vanishes in an asymptotic sense. 



III. COARSENING ERA 
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The largest cluster size N{t) satisfies ODE p.ip . In the 
new units this ODE reads 

N^rjN^-l. (3.1) 

The scaled PDE (|1.3p and conservation identity (|1.6p are 
now 

dtr + dn{{vn^ -l)r^ ^0, (3.2) 



The apparent "road-block" to further growth is not 
the end of the aggregation story. The growth era asymp- 
totics are not uniformly valid as t/N^^^ s-oo. As rj de- 
creases, the exponentially small component s in the full 
convection velocity u ^ rjn^ — s in (jl.3[) gains influence 
until it balances the (now small) ryn a . The critical size 
ri, = (s/ri) , where u — 0, "catches up" with the average 
cluster size and is now near Nq. Clusters smaller than the 
critical size n* shrink, shedding monomers and fueling 
the continued growth of the clusters larger than n* . The 
classic process called coarsening has begun. The charac- 
teristic time of coarsening, to be determined shortly, is 
exponentially longer than the characteristic time [t] N^^^ 
of the growth era. During coarsening, the distribution 
widens and eventually fills the whole range of cluster sizes 
from the (growing) maximal size down to zero. The tail 
of the coarsening era is characterized by convergence to 
one of the self-similar distributions predicted by Lifshitz 
and Slyozov. 



A. Coarsening era scaling 

Relative scaling units^ of time t and super-saturation 
T] follow from the balance of all three terms in the convec- 
tion velocity u = rjn^ — s in ()1.3|) . The balance between 
u and s yields Nq/s as the relative unit of time, while 

balancing rjn^/^ and s gives s/N^^'^ as the relative unit 
of rj. The relative unit R/Nq of r follows from the balance 
of the two terms in the RHS of the conservation identity 
(|1.6p . The relative and absolute units of t, rj, and r 
are summarized in the scaling table: 



in < n < iV, and 
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77=1- 



N 



• dn. 



(3.3) 



B. The determination of the supersaturation 

In the analysis of the creation and growth eras, the 
conservation identity explicitly determines rj from r{n, t). 
In the coarsening era this straightforward approach fails: 

1 /3 

By assuming e ^ 1 we also get s/N^' ^1, hence the 
leading order approximation of (|3.3p is 



N 



nr dn = I. 



(3.4) 



The term containing 77 disappears. Physically, most of 
the available monomers are contained in clusters and the 
super-saturation is vanishingly small. To extract 77 from 
r{n, t) we differentiate p.4p with respect to t: 



N 



Nr{N,t) + / ndtrdn = 0. 



(3.5) 



Next, we substitute N from (|3.ip . and dtr from the con- 
vection PDE (|3.2p into (|3.5p and integrate by parts. After 
some some algebra we find that 77 can be expressed as 



nirdn 



(3.6) 



^ Since the convection PDE 1 11.31 1 and conservation identity 111.61 1 
are non-dimensionalized using nucleation era units [t], [n] from 
HA. II IA.2I 1 for t and n, and e is the unit for rj, dominant bal- 
ances in these equations provide scaling units relative to those 
of the nucleation era. For instance the characteristic cluster size 
relative to [n] is A^o in 112.31 1. and the actual unit of n is No[n]. 



In summary, r{n, t) in < n < N satisfies the intcgro- 
differential equation, consisting of the convection PDE 
p.2p with 77 as in p.6p and N as in (|3.ip . An effective 
initial condition is determined by asymptotic matching 
with the tail of the growth era. At 72 = the convection 
velocity is negative, so a boundary condition there is not 
required. 
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C. Changing variables 



D. Initial conditions and early widening 



The growth of the largest cluster size N{t) with time 
implies that PDE (|3.2p has to be solved on a growing 
interval of n. We simplify the numerical solution by using 
the following change of variables first: 



n 



q{x, t) = Nr{Nx,t). 



(3.7) 



The idea should be clear; The normalized cluster size x 
ranges in the fixed interval (0, 1) and q is the distribu- 
tion of cluster sizes in x— space. Wc multiply r by iV so 
that qdx = rdn. The convection PDE (|3.2p for r(?i, t) 
transforms into an convection PDE for q{x, t), 



dtq + {wq} = 0, 



(3.8) 



in < X < 1. Here, w is the convection velocity in x 

space, 



w 



l(ryiV^(x^-x) + (.T-l)). (3.9) 



Boundary conditions are not required, since w vanishes at 
X = 1 and is negative at a; = 0. Equation (|3.6[) translates 
into a functional dependence of 77 upon N and moments 
of q, 



Io<ldx 
x^qdx 



(3.10) 



The largest cluster size N is easily determined from the 
conservation identity (|3.4p written in terms of q and x: 



1 

N 



X q dx. 



(3.11) 



In summary, both 77 and N are found explicitly from 
q{-, t) on the interval (0, 1), and this makes (|3.8p an ex- 
plicit integro-differential evolution equation for q. It is 
convenient to introduce the moments of q(x, t) (them- 
selves functions of time) : 



Mo = I qdx, Mi= [ x^qdx, Mi 
Jq ^ Jo 

Then ((XTU)) and (jXTTI) become 



^ Ml' Ml' 
3 



X q dx. 
(3.12) 

(3.13) 



and the convection velocity w can be written as 



w 



1) 



(3.14) 



The t limit of the coarsening solution for r(n, t) 
should match distribution (|2.5|) . which characterizes the 
tail of the growth era. Hence, we have the effective initial 
condition 



q{x, 0)=NI](N^ (l^x) 



(3.15) 



The RHS in (|2.5[) is translated into x, q variables. Since 
A^o ^ 1, this initial distribution is a tall spike of height 
A^o^"^ concentrated in a narrow interval of x— values near 
x^l: 0<l-x< 0(iV(7^/^). The initial condition for 
the largest cluster size N (in coarsening units) is N{Q) = 
1. 

To our knowledge, the integro-differential evolution 
equation for q(x, t) does not admit an analytic solution, 
so a numerical solution is sought. From a numerical 
point of view, the tall, narrow initial condition (|3.15|) 
is not desirable for two reasons: First, it is narrow, with 
width proportional to ea, and thus resolving it numeri- 
cally would be difficult (for e ^ 1). Second, this initial 
condition depends on e via the dependence on Nq, thus 
for every e we would need to run the computation again. 
Some preliminary asymptotics fixes both issues and sup- 
plies us with a global solution: As long as the distribution 
remains a narrow spike near a; = 1, and thus the three 
moments — Mg, Mi , and Mi — are all near 1, the convec- 
tion velocity w in (|3.14p can be approximated by 



X3 



1 = ^{x-l) + 0{x-lf 



(3.16) 



near x = 1. The convection PDE (|3.8p with w replaced 
by its linearization (|3.16p can be solved analytically: The 
"early" evolution of q{x, t) based upon the linearized con- 
vection velocity (j3.16p is given by the widening distribu- 
tion: 



q{x, t) = 7Vo^e-*/3j (iv|e-*/3 (1 - x)) 



(3.17) 



This asymptotic distribution matches the effective initial 

condition p.lSp for t = 0, and remains valid as long as 

_ 2 

the "a;— width" remains small, -^e*" ^ 1. 



E. Time-shift and the numerical solution 

The strategy now is as follows: First, we assume that 

our numerical PDE solver accurately resolves a distri- 

2/3 

bution of width 6 with Nq ' <^ 6 1. From the 
e— dependent initial condition (|3.15p . we evolve q{x, t) 
according to the asymptotic solution (|3.17p until the 
"x— width" TVq ^/^e*/'^ achieves the value 5. This hap- 
pens at time 



< = 21ogiVo + 31ogi. 



(3.18) 
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The numerical solver takes over for times greater than 
t in p.lSp . The width S is chosen so that it is much 
larger than the numerical discretization of x, so that the 
solution can be resolved, yet much smaller than 1 so that 
the analytic solution remains valid. 

It is convenient to absorb the e— dependent component 
21ogA''o in (|3.18[) by shifting the origin of time. The 
shifted time is 

t'=t-21ogAfo (3.19) 

and the numerical solver is turned on at shifted time 
t' = 3 log I , with the effective initial condition 

q{x,t') = ^j{^{l-x)), (3.20) 

in < X < 1. As desired, the time-shift produces an 
e— independent initial condition for the numerical solver, 
and thus an e— independent numerical solution. For a 
wide range oi S's in Nq <C (5 ^ 1, the numerical solu- 
tion at fixed t' should be close to the asymptotic solution. 
We use this later (see Fig. [S|) to convince ourselves of the 
numerical solver's acceptable performance. The details of 
the numerical solution are spelled out in the Section fVl 
After finding q{x,t') numerically, we reconstruct 
r(n, t) using (|3.7p : 

r(n,0 = ^9(f,i-21ogiVo). (3.21) 

For t < 3 log 5 + 2 log A^o , we use the asymptotic expres- 
sion (|3.17p for g, and for t > 3 log 5-1-2 log Nq, we use the 
numerical solution. Figure [3] shows the numerical solu- 
tion for the coarsening era at various values of the shifted 
time t' . 
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FIG. 4: The numerical solution at various times as found 
using clawpack. Displayed are snapshots from t' = —10 to 
t' = 2. The solution continues to evolve after t' — 2, as it 
converges to the similarity solution. Prior to t' = —10, the 
solution is described by the (analytic) asymptotic solution. 



The coarsening era solution exhibits three phases: 
widening, transition, and similarity solution (also called 
late stage coarsening). During the initial widening, the 
support of the distribution has not yet reached a; = 0, and 
the fraction of clusters which have dissolved completely 
is negligible (See Fig.[S]). The widening is accurately de- 
scribed by the asymptotic solution (|3.17p . During the 
transition, the support of q reaches down to a; = and 
the smaller clusters start dissolving, so the total density 
of clusters decreases. To resolve this part of the solution 
the numerical solver is required. The solution is shown 
in Fig. m During the "tail" of coarsening we can observe 
the convergence of the distribution to a specific similarity 
solutions of the LS model. 

The three phases of the coarsening era can be seen 
in Fig. [5] which shows the (normalized) distance-^ be- 
tween the numerical solution and the asymptotic solu- 
tion (the solid line), and between the numerical solu- 
tion and the discontinuous similarity solution (the dashed 
line). Initially, the numerical solution agrees with the 
asymptotic solution and the normalized distance is neg- 
ligible. Afterward, in the transition, non-linear effects 
and the non-zero width of the distribution cause a widen- 
ing "rift" between the numerical solution and asymptotic 
one. These non-linear effects also drive the numerical so- 
lution towards the similarity solution (which is described 
in greater detail below), until eventually, the numerical 
solution is almost indistinguishable from one of the sim- 
ilarity solutions. 




-20 -10 10 20 

t' 

FIG. 5: The normalized distance between the numerical so- 
lution and the asymptotic solution given by 13.171 (solid) , and 
the distance between the numerical solution and the discon- 
tinuous similarity solution given by (|4.14p (dashed). 



^ We use , ?i ^ to measure the distance between a nu- 

Jo \a{x)\dx 

mcrical solution n{x) and an asymptotic solution a{x). The nor- 
malization is used because the similarity solution decays to 
as i ^ oo and thus a simple norm might give an impression of 
convergence when there is none. 
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FIG. 6: The three moments Mq, Mi, and Mi calculated for 
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the numerical solution. Around t' = —10, their distance from 
1 is noticeable and the asymptotic solution is no longer valid. 
Around t' — —5, Mo departs from 1 clusters start dissolved 
at X = 0. 



Figure [7] shows P{x) for different values of fi. 



Pl.5- 




IV. SIMILARITY SOLUTIONS 

We briefly review the two-parameter family of LS sim- 
ilarity solutions. 

The convection PDE p.8H3.11"|) admits a separation of 
variables solution: 



(4.1) 



q{x,t')^c{t')P{x). 



We start with the temporal part c(<'): In the ODE ((3?T|) 
for N{t), substitute N = 1/Mi and N^/^rj = M0/M1/3 
as follows from the two equations in (|3.13p . We get 



Ml = Aff • 1 



Ml 



(4.2) 



and equations ((3TT21 [3?T3l [411 lOj) imply an ODE for 
cit'): 



where F and /x are time independent constants defined 
by 



F 



xP dx. 



xi P dx 



The solution of ODE (HH) is 



cit') 



F- (/i- l)(i'-t,)- 



(4.4) 



(4.5) 



Here, tg is a time-shift related to the onset of coarsening. 
It is determined later in the paper using the numerical 
solution and knowledge of the similarity solution. 

Given c(t'), we find the spatial part of the similarity 
solution, P{x). Substituting (|4.ip into the convection 
PDE dSIS]), and using ODE gj]) for c we find an ODE 
for P: 



P^ 
P 



/i(x3 — x) + {x — 1) 



(4.6) 



FIG. 7: Profiles P{x) for various values of /i, 1 < /i < |. The 
profiles are normalized so that P{x)dx = 1. Specifically, 
the values of jj, in the figure are 1.05 through 1.5 in steps of 
0.05. The corresponding orders of contact vary from —0.833 
through (bold) and on to 00 (dashed). A function with order 
of contact p < —1 is non-integrable, and therefore unphysical. 

The parameter fi is related to the order of contact of 
P{x) (with zero) at x — 1. The order of contact is the 
power p so that 

P{x) b{l- xf as X ^ 1", 

for some constant 6 > 0. The super-script (— ) indicates 
that the limit is from below. One sees that 



p ~ lim -^(x — 1) 



(4.3) Substituting into (jiTT]) gives 



5/i — 6 . , , 3p + 6 

p = , or equivalently \i = 



2 



2p + 5 



(4.7) 



(4.8) 



Since the convection velocity w in (j3.14p is regular at 
X = 1, the order of contact of ^(a;, t') at a; = 1 is constant, 
independent of time Q . The coarsening era solution is 
discontinuous at a; = 1, so p = 0, and then (|4.8p implies 
/i = |. We therefore expect the numerical solution to 
converge to the = | similarity solution as t ^ 00 . This 
convergence is verified numerically (see Fig. [S)) 
For /i = |. the similarity solution profile P(x) is: 



P = 



125exp(-y^(coth-XV2T) - tanh-^(2£l^ 



(5 _ 2.2/3 _ 3.1/3^ 



Here, /'(x) is normalized so that 

Pdx = 1. 



(4.9) 
(4.10) 
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FIG. 8: Finding the time-shift ts. From around t' = and on 
the value of t' — .r^,,,. stabilizes on —8.9037. 



We estimate F in (|4.4p for this P{x), numerically. It is 

F w 0.632573 (4.11) 

The dark line in Fig. [7] shows P with /i = |. The dashed 
line shows the Coo solution with /i = |. The only pa- 
rameter that is "matched" to the numerical data is the 
time-shift, tg, found via a simple method described next. 



A. Asymptotic Matching with Coarsening Era 

Finally, we determine the additive time constant in 
(j4.5p by examining the long-time limit of the coarsening 
era solution. By substituting q{x, t') ~ c{t')P(x) with 
c{t') as in ((T^ into (PTT^ for Ah , and setting /i = | , we 
find 



or equivalently. 



Mi(i') = 



ts = t' 



{f - ts) 



M{t') ■ 



(4.12) 



(4.13) 



In order to estimate ts, we calculate Mi{t') from the 
numerical solution and plot the RHS of (|4.13p vs. t' in 
Fig. [5] The horizontal asymptote as t' cx) shows that 
ts ~ —8.9037. In Fig. O we see that the numerical so- 
lution indeed converges to the similarity solution with 
/i = I and ts = —8.9037. From i' = 5 and on the numer- 
ical and similarity solutions are practically indistinguish- 
able. Thus, equations (|4.5i 14.91 14. lip and the value of ts 
imply that the coarsening era solution, q{x, t'), asymp- 
totes to 



625 exp(-2y'| (coth^^ {^21} - tanh"^ 



/21 



0.632573 • {f + 8.9037) (5 - x2/3 - x'^/^Y 



(4.14) 



FIG. 9: The numerical solution q{x, t'), scaled so that Mi — 1, 
at different times. Starting on the right at t' = —10 in a 
narrow distribution, and very close to the similarity solution 
with A* = I (Dark hne) at t' = 4. 



V. METHODS 

Here we describe the method in which we solved the 
non-linear convection PDE for the coarsening era. We 
solve ([XSHXTU)) with the initial condition ([X^ us- 
ing LeVeque's conservation law numerical solver pack- 
age, clawpack @, using the Riemann problem solver 
rpladecon (a Riemann solver for conservative convec- 
tion). Since we expect the solution to start from a narrow 
and tall distribution near x = 1 and widen as t increases, 
we use a non-uniform grid that becomes more dense to- 
wards X = I. Specifically, 



lln 



lOn' 



(5.1) 



where m = 2000 is the number of grid-cells. Notice that 
xq = and Xm = 1- This non- uniform grid is chosen 
so that it has a greater resolution where we expect to 
find the biggest gradients, i.e. near x = 1. The non- 
uniform grid was implemented using a variable "capac- 
ity" in the numerical solver. The capacity of a cell de- 
notes the change in mean value which results from a unit 
flux into the cell. A non-uniform grid can be implemented 
on a uniform grid by giving the computational cells that 
correspond to smaller physical cells a smaller capacity, 
and the opposite for the larger cells, and adjusting the 
convection velocity as needed. For more information on 
implementing non-uniform grid-size in clawpack please 
see 0, Section 6.17]. 

We start the numerical solution at t' = —20 = 3 log 5, 
soS ^ e-20/3 ^ 1.2726 X 10-3. The integrals in ([XT^ are 
calculated as Riemann sums"^ By comparing the results 



as t 



oo. 



We use Riemann sums and not the trapezoidal rule. In clawpack, 
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to those obtained from a finer mesh and smaller 6, we 
estimate the relative error to be ~ 1%. 



VI. CONCLUSION 

The original LS theory with its scale invariance de- 
scribes late stage self-similar coarsening, but not the ac- 
tual process of how it arises from the initial condition of 
pure monomer. In particular, there is no prediction of 
characteristic cluster size, nor characteristic time that 
marks the onset of coarsening. This paper, together 
with its predecessor on the creation era fills this gap in 
a limited sense: It tells the story of the intermediate 
processes according to a "classical" aggregation kinetics 
based on a conservative union of fundamental ideas due 
to Becker-Doring, Zeldovich. and Lifshitz-Slyozov. As 
we have seen, there is a succession of three eras, cre- 
ation^ growth^ and coarsening^ with each consecutive era 
linked to the previous one by asymptotic matching. In 
this paper, there is the additional twist of connecting an- 
alytic solutions to one universal numerical solution for 
coarsening. 

The characteristic cluster size at the onset of coarsen- 
ing is 

Wc = ^N, (6.1) 

where [n] is the characteristic cluster size (|A.2p during 

3 

creation, so the prefactor is the relative increase of 
size between creation and the onset of coarsening. The 
creation era size [n] itself is dominated by the large expo- 
nential expdG^/fc^r), so the free energy barrier against 
nucleation (5* is the most important physical parameter 
in the determination of the characteristic size of ric. 

The time to onset of coarsening is a little more subtle 
than simple scaling alone. From p. 131 14.12|) the largest 
cluster size (in units of [n]c) is 

N{t')^l{t' ~ts)tort'->l. (6.2) 

The long time similarity solution in section lTTTl can be said 
to begin around <' = 0. This is when the graph oft'—jj 
in Fig.[5]becomcs constant, as it should for the similarity 
solution, and when the distance between the numerical 
solution and the similarity solution is small and shrinking 
as can be seen in Fig. [5] Recall that t' is the shifted time 
as in (|3.19p . so equation (|6.2p reproduces the well-known 
result that the largest clusters size grow linearly with 
time during late-stage coarsening. The shifted time t' 
is related to physical time t by translation and scaling, 
according to 

Wc(t' + 21og^). (6.3) 



the cell values represent cell averages, and thus Ricmann sum is 
the correct estimate. 



Here, [t\c is the scaling unit of time from the coarsening 
era, 

Mc = fM (6.4) 
or, using s as in (jA.4p and [t] as in (jA.ip . 

M. = fiw(^). (6.5) 

The factor fs carries the unit of time. Notice that 

[t]c, like [n]c, is proportional to exp(|G*/feBr). This is 
exponentially longer than the creation timescale [t] which 
is proportional to expdC^/fcsT). The additive time 

3 

constant 2 log in (|6.3p comes from (|3.19p with Nq = 

. Setting t' = in (|6.3|) we recover the onset of late 
stage coarsening in physical time, 

tonsot = 2[t]el0g^. (6.6) 

As £ ^ the logarithm term is nominally large. But 
in the "big picture" , it is just another prefactor which is 
overwhelmed by the large exponential exp(| G^/fcsT) in 
[t]c. In summary, the characteristic cluster size and char- 
acteristic time associated with the onset of coarsening are 
both proportional to exp(|G*//sBT). 

VII. DISCUSSION 

Although we have described the three eras, as we have 
set out to do, the aggregation story is not over. According 
to conventional wisdom, the "correct" similarity solution 
is the smooth one with p = oo (and /i = §) and not the 
discontinuous one we found with p = 0. In the original 
LS paper, it is suggested that the (rare) coagulation of 
large clusters eventually selects the smooth similarity so- 
lution. A recent work by Niethammer and Velasquez [3| 
suggests that screening-induced fluctuations also leads to 
the selection of the smooth solution. 

Even classical kinetics without any additional physics 
has structure that is not yet fully examined. It can 
lead to the selection of the smooth solution given suf- 
ficiently long time, much longer than in (|6.4|) . In 
particular, the convection PDE boundary value problem 
for r(n, t) is the lowest order of approximation to a dis- 
crete system of ODE's. The next order of approximation 
introduces an effective diffusion in size space, and this 
smooths out the discontinuity. Another effect of the dis- 
crete kinetics is that the Zeldovich nucleation rate does 
not abruptly "turn on" at t = as we have assumed in 
our reduced analysis. It has been shown (both analyti- 
cally and experimentally [ij ) that there is a transient 
during which the nucleation rate smoothly increases, with 
the Zeldovich rate as its long time limit. This is the so- 
called "ignition transient" . Taking this transient into 
consideration would add a narrow boundary layer to the 
sharp front, connecting it smoothly to zero. A prelimi- 
nary analysis shows that a narrow boundary layer about 
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n = N widens on a timescale much longer than [i]c, and 
the smooth solution results. 

In summary, the discontinuous similarity solution is 
structurally unstable due to a variety of physical and 
mathematical perturbations. It is now a question of time 
scales: The mechanism that causes the fastest devia- 
tion from the discontinuous solution will determine the 
timescale of this last era, and will be the main cause for 
the smoothing of the distribution. This "contest" will be 
played out in a future paper. 



APPENDIX A: COARSENING ERA SCALING 

The scales [t], [r], and [n] of time, cluster density and 
cluster size of the creation era are found (in [l|) to be 

[t] = (87r)-s {eia-i}e^^ {D\f^uj^y^ , (A.l) 



[n] ^ (7rTO2TO\/3) 



(A.2) 



[r] = (3-2"7r^)-i/5 



(A.3) 



Here, e denotes the initial (small) value of the super- 
saturation 77(0). In the exponents, the fraction ap- 
proximates the initial free energy barrier against nuclc- 
ation in units of fc^T. In the prefactors, w is the evap- 
oration rate so that ujns is the rate at which monomers 
at the surface of an n— cluster leave it. The dominant 
balances leading to these scaling units are based phys- 
ically upon the Zeldovich rate of nucleation, diffusion 
limited growth of created clusters, and conservation of 
monomers. In particular the exponential largeness of 
characteristic time and cluster size [t] and [n] in e arise 
from the exponential smallncss of the Zeldovich nucle- 
ation rate (proportional to cxp — (T'^/2?7'^). The relation 
[n] (X [t]2 as evident from (jA.H IA.2P is a signature of 
diffusion limited growth (in 3D). 

The scaled surface tension s in the scaled version of 
the LS equation (|1.5[) is given by 



{i{A^fY {Dv'^ fs) 



M, 



(A.4) 



From (|A.1[ IA.2|) we see that s (x exp(— ■io'^/2£^) is ex- 
ponentially small for e ^ 1. 



APPENDIX B: GROWTH ERA SOLUTION 

Let us, for the moment, take ij{t) as given. The char- 
acteristic curve corresponding to the "first" cluster — the 
one that nucleated at time t ~ — is n = N(t), where 
N{t) satisfies 

N = 7]N^ , 7V(0) = 0, N{t) > for t > 0. (B.l) 



The support of r{n, t) lies in TZ: 

n = {{n,t):0<n<N{t),t>0}. (B.2) 

In TZ the value of r(n, t) is foimd from 

r(n, t) = n^igir), (B.3) 

where ^(t) is the constant value of n^r{n, t) along the 
characteristic curve that has ^(t) = 0. For any point in 
TZ, there is one characteristic curve that passes through 
it, so T in (jB.3P is a function of n and t. Given t = 
T(n, t), (|B.3P is the growth era solution for r(n, t) in TZ. 

The asymptotic determinations of g(r) and r(n, t) are 
simple. It follows from (H?! [BTI) that 



Subtracting these equations gives 

pT(n, t) 



Tj{t')dt'. 



(B.4) 
(B.5) 



r](t')dt'. 



(B.6) 



Characteristics with r = 0(1) are launched during the 
creation era, and for these we have that (?(t) in (|B.3p is 
in fact j(r) from p.9p . During creation, r](t) (in units of 
e) differs from 1 by 0(£^), so for r = 0(1), we replace 
r]{t') in (|R6| by 1, 



(B.7) 



In the limit > 1, the RHS of ((BJ| remains 0(1) for 
N — n = 0{N3) and in this case we replace the RHS by 
its linearization about n = N, so 



N~n 
1— ■ 



(B., 



Once we determine N = N{t), (|B.8P gives r(n, t) and 
the solution for r{n, t) in the "front" 



is given by 



Q < N{t)~n^Q{N{t)Y . (B.9) 



r(n, t)-A^-3jf^V (B.IO) 



For N — n ^ Na, t S> 1 and j(t) asymptotes to zero, 
corresponding to negligible production of new clusters 
after the creation era. 

We complete the story of the growth era by an asymp- 
totic determination of N{t). Since the support of r{Ti, t) 
is effectively the narrow front (|B.9[) . the conservation 
identity p.6p reduces asymptotically to 



1 



N. 



(B.ll) 
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Here, R is the (scaled) total number of clusters produced 
during nucleation, given by p.lOp . Combining (|B.1[ 
IBJT|) gives a simple ODE for N{t), 

iV~(l-f)iV^ (B.12) 

where 

No EE ^. (B.13) 



The solution with iV(0) = (and TV > for t > 0) is 
given implicitly by 

^ = ^r,log(l + r,(f)^). (B.14) 

j=o ^ ^ 

Here, are the cube roots of —1: ro = e* a' , ri = — 1, r2 = 
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